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PREFACE 

The general objective of this project was to test the feasibility of using multi - 
spectral techniques for the determination of water depths from spacecraft altitudes. 
Previous research with multispectral aircraft sensors sponsored by the National 
Oceanic and Atmospheric Agency’s Spacecraft Oceanography Group (formerly 
of the Naval Oceanographic Office), under the direction of John Sherman III, had 
demonstrated that narrow spectral bands in the blue green portion of the spectrum 
could be used for remote bathymetry. 

This report is submitted in fulfillment of NASA Contract No. NASS -2 1783, 

Task I. The principal investigator on the project was Fabian C. Polcyn, and the 
work was carried out by the Infrared and Optics Division of the Environmental 
Research Institute of Michigan under the direction of Mr. R. R. Legault. 

The authors wish to acknowledge the contributions of W. L. Brown who 
worked on the early formulation of the technique and Michael Gordon for 
assistance in developing the computer programs for this task. 
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REMOTE BATHYMETRY AND SHOAL DETECTION WITH ERTS 

1 

INTRODUCTION 

The successful launch and operation of ERTS-1 has provided a wealth of information in 
the form of multispectral imagery and data tapes covering all the continental U. S. and its 
coastal waters, as well as a large fraction of the shallow-water areas in the rest of the world. 

In these shallow-water areas, bottom features can frequently be seen on ERTS imagery to 
depths of nearly 10 meters in clear Caribbean waters and of 2-3 meters in more turbid northern 
waters. In view of the pressing need for updating hydrographic charts in coastal areas, the 
present study was initiated in an effort to extract quantitative water-depth information from 
ERTS data. 

Aerial photography has long been used by the U, S. Coast and Geodetic Survey for the 
purpose of locating aids to navigation, mapping shoreLines, and revealing bottom features to be 
investigated [1, 2]. The U. S, Naval Oceanographic Office also has used ste reo-photogrammetric 
techniques for determining water depths in certain areas [3]. Numerical techniques for calcu- 
lating water depth on the basis of the intensity of bottom-reflected solar radiation have been 
developed at ERIM and successfully applied to aircraft multispectral data [4]. These tech- 
niques now have been applied to ERTS data. In addition, new methods have been developed 
specifically for ERTS data. 

In this report, these techniques are documented and the results of their application to two 
test sites discussed. The use of ERTS data for the purpose of remote bathymetry is then 
evaluated in terms of relative costs and accuracies, as compared with conventional ship-survey 
techniques. Finally, the results obtained from ERTS data are compared with those which might 
be obtained with space sensors that are more optimally designed for this application. 

Two computer algorithms have been developed for extraction water-depth information 
from ERTS data. The ratio method (see Appendix A) corrects for differences in bottom re- 
flectance and/or water attenuation and gives accurate results for depths up to the attenuation 
depth in ERTS Band 5 (about 3 meters in clear water). The optimum decision-boundary method 
(see Appendix B) extracts the maximum amount of depth information in the case of uniform 
bottom reflectance and water attenuation, but it is liable to errors when these parameters are 
not uniform throughout the scene. 

Accuracy could be improved, and changes in bottom topography could be detected, by a 
time study of a sequence of ERTS frames for a given area. It is recommended that a computer 
method of comparing these frames be developed in conjunction with the MIDAS system currently 
being developed at ERIM. 
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In the meantime, processing of ERTS data should continue with presently developed tech- 
niques in areas where doubtful shoals exist or where erosion potential is high. Applications of 
these techniques could contribute interim data valuable to the shipping industry and establish 
baseline information for future erosion studies. 


2 

APPLICATIONS AREA DISCUSSION 
2.1 UPDATING NAVIGATION CHARTS 

The history of navigation and shipping is filled with examples of ship losses resulting 
from collisions with uncharted shoals or shoals whose positions were known only approximately. 
The International Hydrographic Office has expressed concern over the status of navigation 
charts around the world and has called upon the world community to assist in providing more 
reliable chart information. 

One of the benefits that can be predicted from the use of satellites in mapping the oceans 
is improving nautical charts where better location and depth information of shoals would help 
reduce ship casualties. Additional benefits might be in helping to shorten commercial shipping 
routes and to increase safety. The reduction in commercial-ship casualties would reduce 
losses of life, cargo, and property. Table 1 [5] gives an estimate of the average losses for 
the U.S., which run over $90 million per year. Since U.S. shipping is about 10% of the total 
for commercial ships of 100 or more gross tonnage, about $1 billion is the estimated world- 
wide loss. The number of vessels involved and the primary cause of U.S. ship casualties are 
given in Table 2 [5] . At least one loss category is directly attributed to inaccurate depth in- 
formation. In addition, it is possible that other loss categories might be indirectly reduced 
with better navigation charts, which could reduce calculated risks or shorten shipping routes. 

In any event, the reduction in loss on a world-wide basis is estimated at a few million dollars 
per year. 

The use of ERTS remote bathymetry for this purpose has several advantages over conven- 
tional ship-survey techniques. First, the coverage is continuous, providing a two-dimensional 
image rather than a series of points or transects. This means that all shoals are much more 
likely to be detected, if they are within the range of the technique and their geographical lo- 
cations can be easily determined by reference to land areas within the image. Second, ERTS 
coverage is repetitive, covering the same area every 18 days, so that changes resulting from 
coral growth and shifting sand bars can be assessed. Finally, because remote bathymetry is 
less expensive per unit area than conventional methods, areas which were previously uncharted 
because of economic constraints can now be charted. 
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TABLE 1. ESTIMATED U.S. LOSSES THROUGH U.S. COMMERCIAL SHIP 
CASUALTIES (IN DOLLARS) 


Kind of Loss 

1966 

1967 

1968 

1969 

1970 

1971 

Vessel 

95,139,000 

53,080,000 

63,206 >000 

68,267,000 

69,274,000 

78,961,000 

Cargo 

7,454,000 

9,801,000 

5,186,000 

10,269,000 

17,360,000 

6,629,000 

Property 

3,131,000 

12,262,000 

12,676,000 

7,926,000 

10,629,000 

8,911,000 


This table has been constructed using fiscal-year figures published in the following issues 
of the Proceedings of the Marine Safety Council: December 1966; November 1968; December 
1969; December 19J0; and December 1971, 


TABLE 2. NUMBER OF VESSELS INVOLVED IN U. S. COMMERCIAL SHIP 
CASUALTIES — PRIMARY CAUSES 


Primary Cause 

1966 

1967 

1968 

1969 

1970 

1971 

Personnel Fault 

550 

551 

984 

1,193 

1,300 

1,325 

Calculated Risk* 

257 

386 

307 

124 

22 

10 

Restricted Maneuvering Room 

210 

177 

50 

32 

31 

22 

Storros-Adverse Weather 

374 

202 

446 

253 

274 

370 

Unusual Currents 

50 

34 

70 

46 

57 

19 

Sheer, Suction, Bank Cushion 

36 

31 

34 

57 

3 

4 

Depth of Water Less than Expected 

100 

63 

110 

76 

54 

50 

Failure of Equipment 

311 

376 

540 

610 

510 

524 

Unseaworthy-Lack of Maintenance 

340 

323 

135 

133 

86 

81 

Floating Debris-Submerged Object 

153 

157 

97 

151 

172 

151 

Inadequate Tug Assistance 

137 

105 

34 

63 

30 

14 

Fault on Part of Other Vessel or Person 

746 

947 

1,142 

1,333 

1,304 

1,435 

Unknown -Insufficient Information 

29 

23 

62 

112 

220 

147 

TOTAL 

3,293 

3,030 

4,001 

4,183 

4,063 

4,152 


The source of the figures is the same as for Table 1, 

*Before 1969, "Error in Judgment-Calculated Risk," rather than just "Calculated Risk.” 


2pt 


FORMERLY" WILLOW RUN LABORATORIES, THE UNIVERSITY OF MICHIGAN 


In order to take best advantage of the greater range and accuracy of conventional tech- 
niques, the two methods should be used in conjunction; satellite bathymetry can be used to 
detect and locate shoals, and conventional ship surveys to confirm their depths at sampled 
points and chart surrounding areas to greater depths. 

2.2 MONITORING SHORELINE CHANGES 

The bottom topography of coastal areas is constantly changing as a result of the natural 
processes of erosion and deposition. These processes are accelerated by storms and changing 
water levels in inland lakes. They are also modified, either intentionally or unintentionally, by 
the construction of harbors, recreational beaches, power-generation plants, marinas, and 
local erosion-abatement devices. 

There is need for a rapid monitoring technique to measure changes in coastal topography 
over wide areas. Such a capability could be used to assess storm damages, identify areas of 
high erosion risk, and evaluate the impact of shoreline construction. Information gathered by 
this process could be used, for example, by the Environmental Protection Agency for zoning 
purposes or by the U.S. Army Corps of Engineers for building more effective erosion-control 
devices. 

Because of its repetitive wide-area coverage and low cost, ERTS obviously has potential 
for this type of application. Aircraft multispectral data and conventional depth- sounding 
techniques could be used in conjunction with ERTS data to obtain greater accuracy and resolu- 
tion. The techniques of remote bathymetry are documented in the following chapter, and the 
potential of ERTS data for this application is evaluated in the remainder of this report. 


3 

EXTRACTION OF WATER- DEPTH INFORMATION FROM ERTS MSS DATA 
3.1 GENERAL CONSIDERATIONS 

The ERTS multispectral scanner receives electromagnetic radiation in the four spectral 
bands indicated in Table 3. The radiance at the detector is measured and converted into an 
integer data value between 0 and 63. In the linear mode the signal is directly proportional 
to the radiance at the detector. That is, 

Vi = ^ (1) 

where is the radiance in band i at the detector, and is the sensitivity constant for the 
detector (see Table 3). In the signal- compression mode of ERTS, the signal is not linearly 
proportional to the radiance (see Figure 1) but the relationship may be considered approxi- 
mately linear over the middle range of signal values. 
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TABLE 3. ERTS MSS-DETECTOR CHARACTERISTICS [6 j 


Band 

Wavelength Range 
(4m) 

Maximum Radiance 
(mw cm'2 sr‘*) 

Sensitivity Constant 
(mw~l cm^ sr) 

4 

0. 5-0.6 

2.48 

25.4 

5 

0.6-0. 7 

2.00 

31.5 

6 

0,7-0. 8 

1.76 

35.8 

7 

0. 8-1.1 

4.60 

13.7 
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The radiance observed over shallow water is the result of sunlight reflecting from the 
bottom and the water surface, as well as of the scattering of sunlight in the water and the at- 
mosphere. That part of the signal resulting from bottom reflection contains information about 
the depth of the water through which the light has passed. In order to extract this information, 
one must first separate the bottom -reflection signal from the rest of the observed signal, and 
then determine how this signal is related to the water depth. 

Assuming the linear relationship (see Eq. 1) is correct, the signal observed over shallow 
water may be expressed as follows: 


V, = k.R . + k.r.R . 
1 i pi 1 1 si 


+ k.r. * 

i i 


H i i - (sec 0 + sec <j> )z 

* V 1 


where R . 

Pi 

T. 

1 

R . 
si 

H. 

l 

p i 


c*. 

i 

Q 

* 

z 


= atmospheric path radiance 

- atmospheric transmittance 

- radiance resulting from (water) surface reflection 
= solar irradiance incident on water surface 

= bottom reflectance 

- index of refraction of water 

- attenuation coefficient of water 

= angle of observation (under water) 

- solar -zenith angle (under water) 

= water depth 


( 2 ) 


This equation neglects volume scattering of sunlight from the water itself, which is usually 
small compared to the other components when conditions are favorable (i.e., when the water 
is fairly clear). 

The first term of this equation accounts for the scattering of sunlight by the atmosphere; 
the second for the specular reflection of diffuse sky radiation by the water surface. The re- 
flection of direct sunlight from the water surface into the scanner, known as sun glint, is 
avoided by restricting angles of observation to less than the solar -zenith angle (see Figure 2). 


Since these two terms contain no information about the water depth, they are considered 
as the background signal which must be removed before the actual depth processing is done. 
This background signal, 


V. . =k.R . + K.r.R . 
bi i pi i i si 


(3) 


may be determined by scanning over deep water, where there is no bottom-reflected signal. 
When this background signal is subtracted from the data, the remainder represents the light 
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FIGURE 1. MSS OUTPUT COUNT (Vi'i VERSUS 

RADIANCE (R./R „ ), COMPRESSED AND 
l max 

LINEAR MODES 



FIGURE 2. GEOMETRY OF WATER-DEPTH MEASUREMENT 
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which has been reflected from the bottom and attenuated in the intervening water layer. From 
Equation (2), this part of the signal may be expressed as 


H. . -a (sec 0+sec <i> )z 
AV i = VjL - V bi = ViT p i ^2 e 

n 


( 4 ) 


The depth dependence of this signal arises from the attenuation of light in the water, expressed 
by the exponential factor in this equation. Methods of calculating the water depth at each point 
in the scene are discussed in the following four sections. 


3.2 SINGLE-CHANNEL METHOD 

If all the parameters in Equation (4) are known, it may be inverted to yield the water depth, 
as follows 


z 


1 

a. (sec & + sec 0 ) 



where 


(5) 


AV 


01 


H. 

k. r. — rP . 


1 1 77 


l n 2 


( 6 ) 


This latter quantity, which is interpreted as the bottom-reflected signal at zero water depth, 
may be determined either from Equation (6) or from examining the scanner signal at the 
water's edge. 

The attenuation coefficient a. may be estimated from Figure 3, although there is a wide 
variation in the attenuation coefficients depending on the purity of the water. Alternately, if 
the water depth is known at one point in the scene, the attenuation coefficients may be deduced 
from Equation (5). They may also be measured in the laboratory, of course, if samples of the 
water are taken from the scene. 


With ERTS data, the scan angle 9 is small enough so that sec 9 may be approximated by 1. 
There is then a unique signal level associated with each water depth, assuming that the water- 
attenuation coefficient and bottom reflectance remain constant throughout the scene. A water- 
depth chart may then be produced simply by mapping the most penetrating channel (i.e., the 
channel with the smallest attenuation — channel 4 for ERTS). In such a map, a distinct symbol 
is printed out at each point in the scene for each range of data values specified. By using 
Equation (5), we can specify these ranges to correspond to convenient depth ranges. 

The advantage of this method is its simplicity and relative insensitivity to random noise 
in the data. It has two drawbacks, however. First, from an operational standpoint, it is 
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difficult to accurately estimate the water -attenuation coefficient without some ground-truth data: 
either the depth must be known at one point in the scene, or a sample of water must be collected 
and analyzed in the laboratory. Second, if there are variations in either the water attenuation 
or the bottom reflectance within the scene, systematic errors will be introduced into the calcu- 
lated depth. Both of these drawbacks are eliminated or reduced — but at the cost of greatly 
increased sensitivity to random noise — in the ratio method discussed in the following section. 


3.3 RATIO METHOD 

If one takes the ratio of the bottom -reflected signals in two spectral bands i and j, Equa- 
tion (4) yields: 


AV i k i r i H i^i -(o^-ttjXsec Q+ sec <j >) z 

AV7 = k.r.H.p. e 
J J 3 J J 


(7) 


This equation may be inverted to give the following equation for water depth: 

1 Ml \ 

z _ * f n j 1 R 1 

(a. - a.)(sec 9 + secT^ Uv. jiJ 


( 8 ) 


where 


k.r.H.p. 
r = 3 } J 3 
ji k i r i H i P i 


(9) 


This method of calculating water depth from two channels has several advantages over the 
single-channel method. As mentioned in the previous section, the attenuation coefficients and 
bottom reflectances change from one location to another, from one time to another at the same 
location, and even from one point to another within the same scene. However, the difference 
between the attenuation coefficients in two appropriately chosen bands exhibits much less vari- 
ation from point to point and from time to time {see Figure 3). The same is true of the ratio 
of the bottom reflectances in the two channels. Since only this difference and this ratio appear 
in Equation {8), the results using this method are much less sensitive to changes in water 
quality and bottom type than are the results of the single-channel method. This means that (1) 
it is easier to specify input parameters for this method, and that (2) there are fewer systematic 
errors inherent in the results. 

The ratio method of extracting water depths from multispectral data appears to be a 
powerful and useful technique. However, it does require optimum channel selection and a high 
signal-to-noise ratio in the data. ERTS bands 4 and 5 are useful for water-depth work but not 
optimum. In addition, the low gain setting reduces the usefulness of channel 5 for determining 
maximum water depths. Even under the best conditions, the bottom -returned signal in channel 5 
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falls below the noise level for depths greater than a few meters* The non-linear response 
characteristics of ERTS in the compressed mode also cause problems in the application of 
this method to ERTS data. 

In processing the Bahama Bank data, a two-step procedure was used to extend the range 
of the depth chart. First, the ratio method was applied, using minimum coastal values for the 
water -attenuation coefficients (from Figure 3). This resulted in a chart with a maximum 
depth of about 2,8 meters. These depths were then visually correlated with the channel- 4 data 
values, and this correlation was then used to determine the actual attenuation coefficient for 
channel 4. By this procedure we determined that the water actually falls in the mean oceanic 
category (the ratio- depth results are still valid for this case because the difference between 
the channel-4 and channel -5 attenuation coefficients is the same for both types of water). This 
attenuation coefficient was then used to produce a single -channel depth chart, with a maximum 
depth of about 9 meters. 


3.4 OPTIMUM DECISION-BOUNDARY METHOD 

In processing the Lake Michigan -Green Bay Inlet data, the ratio method was found to be 
restrictive because of the poor quality of the data. The combination of low sun elevation (30°) 
and relatively high water attenuation restricted the penetration depth to about 2 meters in 
channel 4 and less than 1 meter in channel 5. The number of points for which depth could be 
calculated by the ratio method was too small to apply the two-step procedure used for the 
Bahama Bank data. Consequently, the data were examined more thoroughly, and a third method 
was developed for extracting water depth from this type of noisy data. 

Equation (4) predicts a correlation between the data values observed over water in chan- 
nels i and j, of the form 

a. 

log (AVj) = log (AV.) + log k (10) 

1 J 


where k is a constant depending on bottom reflectance, scene illumination, etc. Points se- 
lected from the Lake Michigan test site are plotted in Figures 4 and 5. The log -log plot 
(Figure 5) shows a linear correlation between AV^ and AV^, with a regression coefficient: 



(ID 


The ratio of the attenuation coefficients in the two channels can be obtained by this method. 
Therefore, if the difference in the attenuation coefficients can be estimated from Figure 3, 
their absolute values can also be determined without ground -truth data. 
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Figure 5 also suggests a method for estimating depth which is less sensitive to noise than 
the single -channel method. In the single-channel method, the decision boundaries for depth are 
vertical lines in Figure 5, However, assuming the actual data points are randomly distributed 
around the orthogonal regression line in Figure 5, the optimum-decision boundaries for classify- 
ing these points are lines drawn perpendicular to this curve. Thus, if a given data point falls 
between the lines marked z = I m and z = 1.5 m, the water depth at this point should be classified 
in this range. This procedure is mathematically implemented by the following equation: 

1 N /AV A 

z -cIJV oe (ivf) (12) 

where 

N 2 

c = (sec 0 + sec 0)^^ oC (13) 

i=l 

and AV^ is defined by Equation (6). 

This method may be applied to any number of channels (N). The main advantage is an 
improvement over the single-channel method in the accuracy of classification* Variations in 
water attenuation and bottom reflectance may cause systematic errors, as in the single- channel 
method, but they can perhaps be minimized by a judicious choice of channels. These errors 
also can be reduced by means of the time-study method described in the following section. 

Equation (12) was implemented in a computer program (see Appendix B) and applied to 
the Lake Michigan test site, using ERTS bands 4 and 5. No attempt has yet been made to find 
the optimum channels for this method, since only two penetrating channels are available on 
ERTS. Results of this application are discussed in Section 4. 

3.5 ERROR REDUCTION BY TIME-STUDY TECHNIQUES 

Errors introduced by variations in water quality, clouds, and other transient phenomena 
can be reduced by a repeated analysis of the same scene over a period of time, using any of 
the three methods described above. Averaging the results on a point- by -point basis should 
reinforce the stable features, as well as reduce the transient effects. 

A photographic overlay method would probably be the most feasible way of carrying out 
this averaging process at present. While the capability of averaging two scenes on a point-by- 
point basis by computer is not yet available, it should be possible in the near future with the^ 
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MIDAS hybrid computer system now under development at ERIM. The repetitive coverage 
afforded by ERTS also should qualify ERTS data for this type of treatment. 


4 

DISCUSSION OF RESULTS 

4.1 LITTLE BAHAMA BANK 

The primary test site for this study had been designated as an area west of Puerto Rico 
where several doubtful shoals are reported on the Hydrographic Office charts. However, 
this area was obscured by clouds during the first ERTS passes, and no useful data were ob- 
tained. Consequently, processing was started on an alternate test site north of Grand Bahama 
Island, covered by ERTS frame 1079-15165 (see Figure 6). 

These data w T ere obtained on 10 October 1972 at 16:15 GMT, when the solar elevation was 
47.4°. The principal point of the image is at 27.45° N, 78.82° W (northwest of Grand Bahama 
Island). The frame covers approximately 100 x 100 nautical miles, with Grand Bahama Island 
in the lower right-hand corner and the Little Bahama Bank covering most of the lower half of 
the frame. There is a 30% cloud cover, but most of the clouds are in the upper half of the 
frame, over deep water. There is a small patch of clouds in the center of the Bank and over 
Grand Bahama Island. 

MSS band 4 (0.5 to 0.6 \i m) shows many underwater features north of Grand Bahama 
Island. Most of these can be identified on the Hydrographic Office depth chart (see Figure 7), 
but a number of streaks radiate from the center of the Bank which are not on the chart. The 
significance of these streaks is discussed below. 

MSS band 5 (0.6 to 0.7 |im) shows some underwater features, but the penetration depth of 
this channel is clearly less than that of MSS band 4. 

MSS band 6 (0.7 to 0.8 /im) and 7 (0.8 to 1.1 fim) show no underwater features, but they 
are useful for defining the land -water interfaces. 

Digital processing of this data began with our applying the ratio method described in 
Section 3.3 to MSS bands 4 and 5. The computer program implementing this method is listed 
and described in Appendix A. The input parameters for this program were obtained from the 
Smithsonian Physical Tables, except for the background or deep-water signals which were ob- 
tained from the scanner data itself. These values are shown in Table 4. 

The maximum depth which can be computed by this method is the depth at which the 
channel- 5 signal falls below the noise level. For this data set, using only the best of the six 
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ERTS detectors, this maximum depth is about 2.8 meters. The depth chart produced by this 
method is shown in Figure 8. 

In order to extend the depth range of this chart, a second map was made using only chan- 
nel 4, which penetrates to a much greater depth than channel 5. To calibrate this single- 
channel map, the channel- 4 data values were correlated with the depths obtained from the 
ratio method. From this correlation, an attenuation coefficient of 0,10 m was inferred for 
channel 4. This is in agreement with published values for mean oceanic conditions (see Fig, 3). 
Using this attenuation coefficient, we extrapolated the relationship between the channel-4 sig- 
nal and the water depth to a depth of 9 meters. This relationship is shown in Table 5, and the 
depth chart based upon this relationship is shown in Figure 9. 

There are a few discrepancies between the depths indicated on these two maps, notably in 
the area around Memory Rock (on the western edge of the Bank) where the ratio method gives 
greater depth values than the single-channel method, and around Cormorant Point (on the cen- 
ter of Grand Bahama Island) where the reverse is true. These discrepancies probably are due 
to the presence of different bottom materials in these areas. The area around Memory Rock 
appears bright on the MSS-4 imagery and dark on the MSS- 5 imagery, indicating that the bottom 
material is probably a bright green coral. The area around Cormorant Point is dark on both 
the MSS-4 and MSS- 5 imagery, perhaps the result of dark organic material washed from the 
land. In both cases the ratio method gives more accurate depths, as expected. 

The single-channel map (Figure 9) shows two other anomalies (see Fig. 9b), First there 
are erroneous depth values near the clouds in the center and at the upper edge of the scene. 
These are due to cloud shadows and scattering from the clouds themselves. Second, there are 
numerous streaks apparently radiating outward from the center of the Island, where the depth 
is indicated between 1 and 2 meters. These could be either sediment plumes formed by schools 
of fish, or shifting sand bars. Aerial photographs of both these phenomena (Figures 10 and 11) 
confirm the former hypothesis. ERTS frames taken on 8 January 1973 (Figure 12) and 26 
April 1973 (Figure 13) show streaks of varying length and direction similar to those in 
Figure 10 but in different locations thus supporting their association with the action of fish 
schools in search of food. 

4 .'2 LAKE MICHIGAN 

The second data set processed during the course of this study was ERTS frame 1089-16090, 
which covers the northwestern corner of Lake Michigan, including the inlet to Green Bay where 
there are numerous shoals and small islands, and the east coast of the Door Peninsula where 
there are several shallow bays (see Figure 14). The date of this ERTS pass was 20 October 
1972, The solar elevation wa&30° during the pass, resulting in a scene illumination almost 
50 percent lower than in the Bahama Bank scene. 
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FIGURE 8. RATIO DEPTH CHART, LITTLE BAHAMA BANK 
Original figure in color. 
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TABLE 4. INPUT PARAMETERS FOR DIGITAL 
PROCESSING OF BAHAMA BANK DATA 

Input 


Parameter 

Value 

<V a l 

0.26 m" 1 

p 2 /p l 

1.26 

V H 1 

0.92 

T 2 /T l 

1.07 

\ 

22 counts 

\ 

11 counts 

0 

42.6° 


TABLE 5, CHANNEL- 4 SIGNALS 
VERSUS DEPTH 


Channel- 4 



Signals 

Depth 

Chart 

(counts) 

(m) 

Symbols (Fig. 9) 

20-24 

9-oo 

W (blue) 

25-26 

7 - 9 

X (blue) 

27-29 

5 - 7 

* (blue) 

30-33 

3 - 5 

• (blue) 

34-36 

2 - 3 

= (red) 

37-40 

1 - 2 

X (red) 

41-45 

0-1 

1® (red) 
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(a) Large -Area Map 

FIGURE 9. SINGLE -CHANNEL DEPTH CHART, LITTLE BAHAMA BANK 
Original figure in color. 
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(b) Insert, showing Anomalies in Depth Values 
FIGURE 9. SINGLE -CHANNEL DEPTH CHART, LITTLE BAHAMA BANK 
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FIGURE 10. 


AERIAL PHOTOGRAPH OF TURBIDITY PLUMES FORMED BY ACTION OF 
FISH SCHOOLS. Little Bahama Bank, November 1973. 
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FIGURE 11. AERIAL PHOTOGRAPH OF SAND BARS, LITTLE BAHAMA BANK. 

November 1973. 
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FIGURE 12. ERTS FRAME 1169-15171-4 (LITTLE BAHAMA BANK, 8 JANUARY 1973) 
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FIGURE 13. ERTS FRAME 1277-15174-4 (LITTLE BAHAMA BANK, 26 APRIL 1973) 
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FIGURE 14, PORTION OF LAKE SURVEY CHART 70, NORTHWEST LAKE MICHIGAN 
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Some underwater features are visible on the MSS band-4 imagery along the Door Peninsula 
and the inlet to Green Bay {Figure 15). Very few underwater features are visible in band 5. The 
water seems to be fairly uniform, but much more turbid than the Bahama -Bank water. 

Because of this combination of low sun angle and high water attenuation, the ratio -processing 
method did not produce any useful results. In an effort to improve upon the single -channel meth- 
od, we examined the data and developed a third processing technique. This method makes 
use of optimum -decision boundaries for depth classification, as described in Section 3.4. From 
the correlation of channel- 4 and channel- 5 signals (see Figure 5) and published water -attenuation 
curves (see Figure 3), the water -attenuation coefficients were taken to be 

a, = 0.30 m * 

1 (14) 

- 0.45 m 1 

Using the combination of channel- 4 and channel-5 signals indicated in Equation (12), the 
minimum signal that can be distinguished from the background noise corresponds to a depth 
of about 2 meters. 

Two depth maps were produced by this method. The first map (Figure 16) covers most of 
the Green Bay inlet, from Point Detour to the Porte Des Morts Passage. The accuracy of this 
map was checked by comparison with Lake Survey Chart 702. Most of the features in the lower 
half of this region are accurately shown in the digital map, but some errors are present in the 
upper half as a result of increased turbidity in the water coming from Big Bay De Noe. Fisher- 
man Shoal, east of Washington Island, is clearly shown on the ERTS map, and the water depths 
around Washington Island agree quite well with the Lake Survey Chart. St. Martin's Island 
Shoals (south of St. Martin Island) are beyond the depth range of the map, as Is Poverty Island 
Shoal (northwest of Poverty Island). Little Summer Island Shoal should be just on the edge of 
detectability, but it is not indicated on the map. A number of false depth values appear in the 
1. 5-2.0 meter range as a result of turbidity to the north of St, Martin Island. 

The second map (Figure 17) covers the eastern coast of the Door Peninsula fi'om Bailey's 
Harbor to the end of the Peninsula. The depths indicated on this map are in good agreement 
with Lake Survey Chart 702. Rowley Bay, North Bay, and Moonlight Bay are accurately mapped, 
and Four Foot Shoal is detectable in the mouth of Rowley Bay. Nine Foot Shoal (near the tip of 
the Peninsula) is beyond the range of detectability and is not shown on the map. 

The depth range of these maps should be at least doubled by using ERTS data with higher 
sun elevation. Data tapes and imagery have been received from ERTS passes over Lake Mich- 
igan during the summer of 1973 which show bottom features to a depth of 4-5 meters. Over a 
period of a few years, repeated ERTS coverage of the Great Lakes should yield data of sufficient 
quality to update navigation charts over large areas. 
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(a) Band 4 

FIGURE 15. ERTS FRAME 1089-16090 (UPPER LAKE MICHIGAN, 

20 OCTOBER 1972) 
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(b) Band 5 

FIGURE 15. ERTS FRAME 1089-16090 (UPPER LAKE MICHIGAN, 
20 OCTOBER 1972) (Concluded) 
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(a) Digital Depth Chart (Strip 1) 

FIGURE 16. COMPARISON BETWEEN DIGITAL DEPTH CHART 
(STRIP 1) MADE BY OPTIMUM DECISION-BOUNDARY METHOD 
FROM ERTS FRAME 1089-16090 AND LAKE SURVEY CHART 702 
(WASHINGTON ISLAND TO POINT DETOUR). The density of com- 
puter symbols denotes depth. Original figure in color. (Continued) 
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OPTIMUM DECISION-BOUNDARY METHOD FROM E RTS FRAME 1089-16090 AND LAKE 
SURVEY CHART 702 (WASHINGTON ISLAND TO POINT DETOUR) (Concluded) 
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(a) Digital Depth Chart (Strip 2) 

FIGURE 17. COMPARISON 
BETWEEN DIGITAL DEPTH CHART (STRIP 2) 
MADE BY OPTIMUM DECISION-BOUNDARY 
METHOD FROM ERTS FRAME 1089-16090 AND 
LAKE SURVEY CHART 702 (MOONLIGHT BAY TO 
ROWLEY BAY). The density of computer symbols 
denotes depth. Original figure in color (Continued) 
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to Rowley Bay) 

FIGURE 17. COMPARISON BETWEEN DIGITAL DEPTH CHART (STRIP 2) MADE BY 
OPTIMUM DECISION-BOUNDARY METHOD FROM ERTS FRAME 1089-16090 AND LAKE 
SURVEY CHART 702 (MOONLIGHT BAY TO ROWLEY BAY) (Concluded) 
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5 

EVALUATION OF EHTS POTENTIAL FOR REMOTE BATHYMETRY 

5.1 RELATIVE COSTS AND ACCURACIES OF REMOTE SENSING VERSUS SHIP SURVEYS 

Costs for processing ERTS data to produce depth charts have been estimated to be on the 

order of $1.50 per square mile. Total costs for collecting and processing multispectral data 
are difficult to estimate, but they are probably at least an order of magnitude smaller than 
those involved in conventional ship-survey techniques. Accuracies of remote-sensing tech- 
niques are on the order of 10-20 percent for depths at which multispectral techniques can be 
used to correct for differences in bottom reflectance and water attenuation (i.e., where bottom- 
reflected signals are observed in at least two channels). 

The accuracy of modern ship- survey techniques is on the order of 2-3 percent. Some 
charts, however, contain data based on survey records from the early 19th century, when the 
available techniques probably yielded less accurate results than present-day, remote-sensing 
techniques. The time and equipment required for conventional ship surveys have prevented 
many countries from updating these charts. Even the maps distributed by U. S. agencies con- 
tain many doubtful shoals indicated by ’’position approximate” and rr existence doubtful.” 

ERTS data is presently available, and is still being collected, for areas in which conven- 
tional ship- survey techniques are too slow and too costly to adequately update hydrographic 
charts or to keep up with changing conditions. In these areas, the use of remote bathymetric 
techniques in conjunction with ERTS data could make a valuable contribution to the store of 
hydrographic information now available, 

5.2 EVALUATION OF ERTS VERSUS BETTER SPACE SENSORS FOR MEASURING WATER 

DEPTH 

The utility of ERTS data for bathymetric purposes is limited by the gain and band location 
of the ERTS sensor. With the low gain setting normally used, an approximate dynamic range 
of only 20 counts can be observed in the bottom -reflected signal under optimum conditions. 

The water-depth range corresponding to this signal range is approximately equal to the atten- 
uation length 1/a. This depth range could be increased by almost 50 percent by increasing the 
gain by a factor of 4 and encoding the data to 8 bits instead of 6. 

For clear Caribbean water, the attenuation length is about 10 meters in ERTS band 4 and 
about 3 meters in band 5. These could be increased by narrowing the bands and positioning 
them closer to the wavelength of minimum water attenuation (c.f., Figure 3), 

An example of a better space sensor for measuring water depth is the Skylab S-192 multi- 
spectral scanner. 
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Band 3 (0.50-0.55 pm) of this sensor has half the bandwidth of ERTS band 4, and it is po- 
sitioned very near to the point of minimum water attenuation. The attenuation length of Car- 
ibbean water in this channel is about 12.5 meters. This, combined with the 8-bit precision of 
Skylab data, allows depths of 15 meters to be seen on Sky lab scanner data. 

The repetitive coverage of ERTS is a positive feature which might be explored further, 
using the time-study techniques described in Section 3.5. 

These could be used to reduce errors and study changes in bottom topography over a period 
of time. 

In summary, although ERTS is not an optimum sensor for this application, there is a large 
amount of depth information in ERTS data which could be extracted by the techniques documented 
in this report. The value of this information would seem to justify further processing of ERTS 
data by these techniques, as well as developing new techniques for comparing a time sequence 
of ERTS frames of a given area. 
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Appendix A 

COMPUTER PROGRAM FOR RATIO METHOD 

A computer program implementing the ratio method discussed in Section 3.2 was written 
during the course of this contract. This program performs the functions of subtracting the 
background signal (see Equation 3) and computing the depth from a specified pair of channels 
(see Equation 8). This computation is done at each point in the scene, and the results are 
written on an output tape. The program also can compute the depth simultaneously from several 
different pairs of channels. The result of the computation for each pair of channels is written 
out as a separate channel on the output tape, followed by the mean and standard deviation of the 
values calculated from each pair of channels. 

The input variables for this program are summarized in Table A-l. These are specified 
in the flexible READ AND PRINT DATA format in MAD language. The actual computer coding, 
in MAD language, is listed in Table A-2. 
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TABLE A-l. VARIABLES FOR DEPTH1 PROGRAM 


Card 1: 

Input Variables Mode 

NV Integer 

NC Integer 


Card 2: 

Process Variables 


Card 3: 


INBIN 

Integer 

OUTBIN 

Integer 

UNIT 

Integer 

OUNIT 

Integer 

NS A 

Integer 

Input Variables 


COEFSW 

Integer 


DEEP (I) Integer 

NPAIRS Integer 

PAIRS (I) Integer 


ALPHA (I) 
RHO (I) 
TAU (I) 

EK (I) 

HU) 


Floating Point 
Floating Point 
Floating Point 
Floating Point 
Floating Point 


SDA 

DEL 

XMID 


Floating Point 
Floating Point 
Integer 


Description 

No. of channels on input tape 
No. of channels on output tape 
(NC-NPAIRS + 2) 


Bin no. of input tape 

Bin no. of output tape 

Unit on which input tape is mounted 

Unit on which output tape is mounted 

Array containing line and point numbers of scene 


If COEFSW“0, a coefficient tape containing the 
deep-water (surface -reflection) signal at 
each point number must be put in 
If COEFSW=l* a single value for the deep-water 
signal is read in for each channel (see below) 

The deep-water signal in channel r (1=1 , . . NV) 

No. of pairs of channels used 

Channel numbers used (1=1 . . . 2*NPAIRS) 

Water attenuation coefficient in channel I 
Bottom reflectance in channel I 
Atmospheric Transmissivity in channel I 
Scanner sensitivity constant in channel I 
Solar ir radiance at surface in channel I 

Solar -zenith angle 

Angular separation of points along scan line 
Point number of nadir 


VF 


Integer 


EDITCH 

LO 

HI 


Integer 

Integer 

Integer 


Amplitude of fluctuations in deep-water signal 

Channel No. used for editing out land 

Minimum signal observed over water in channel EDITCH 

Maximum signal observed over water in channel EDITCH 
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TABLE A-2. LISTING OF DEPTH1 PROGRAM 

C/ORNAL I IHCT [UN L»£pnu. 


* $ 

* WAV?P'< DEPTH SUiJkUUTI Mt PACKAGE — 11 K A ?< C H 1973 — MF G * 

* ^ 

* mis package rs hcsToned to cal cula r e the depth of * 

* SHALLOW *iA TE R ARE&S BASlO UN Tht DIFFERENCE U FT* FEN * 

* I HE RLE Ltd Alices OS< KADIANCrS CVlR S HALLO* and over * 

* L : F E P N A T F P A 7 L A S « * 

* * 




normal miue is integer 

REE L <CNCE S UN 

fcSASABLt SK P2I22B), OUNITp OKEELt GF1LE, UL 1 Nt * T/JR I T E * OCHAN * 

1 PNWDSi PR GO , MODEL* MUUF2* U.M IT* CAL I, Nit* CI'LINC* NS A , NStt, 

2 KS, NA t NH, KP, IO(U* TFLAG, [PACK* RE3E;<V<4), QFACTKI49J, 

3 C T I T L 7 ( 1 0 1 i 3 UTLl(L 9J, C L l S T l l 9 ) , GSPARG<46), OFLAG* 

4 COA.NG i DBA :\ij * wKP* CNA* CUSS, UNChAN, Ci^GDEt C«ECA* QAECD * 

5 OF ILL* URE EL » tiNWflDS 

ERASABLE: U A T A ( 42 3 ) * ITFST , USTAKT* USTAKT* NV, NX, NC* L, IP* 

1 1CP>* NOP* i\CXT* START, REACT . RFAOL* R c 4DP , LABEL * 

2 F X TRA l 2 0 > * DATUM 124), ICDU£124> # 1CHANU4)* CGEF(l) 


** DIMENSION AND MODE DECLARATIONS ** 

DIMENSION DEEP (24) 

DIMENSION AC OFF ( 400*12 I 

U i F ENS [CN I A { 2) , PAi RS 118*2 > tPMUOEl U *KET ( 24 ) ,MAX ( 24 I 

FLOATING 301 NT F ( 1000) * H ( 24 ) * VS S ( 2 4 ) , ALPHA ( 24 ) , TAIJ ( 24 ) , G( 18 )* At IB*2>* 
1 RHGL24) ,FKAX(24> ,RHLSS(24) ,Cf ACTR 

FLOATING POINT OH 4 -NO * tfl) ANG * SUA , DEL * W* X * P i«C * C 5 1 
FLOATING POINT SCAT . , CUS. ■ ELUG. 

FLOATING POINT t K < 2 4 ) 

BOOLEAN BRMCmb 

** VECTOR VALUES ** 


VECTOR VALUES Ek . . . EK U4 > * l .0 
VtCTOR V At Ut S C4 1 =5 ) *0 
VECTOR VALUES P 1 8 0= 0 - 0 L 74^ 329 

VECTOR VALUES REKR = *H ♦ ORE TURN CODE 0F*,I4,H # FROM COEFFICIENT FILE READ 
1 '*IL*t 

WHENEVER NE X T * G • 3 * FUNCTION RETURN 
TRANSFER TO STEP! NEXT > 

** STEP 1 - PROGRAM INITIALIZATION ** 



45 



FORMERLY WILLOW RUN LABORATORIES, THE UNIVERSITY of MICHIGAN 



STEP! 1 ) CHIfiil = -1 

*ERG. 1CUNIT, 601 rCH f H| ,LQ) 

SP'UY,(1 m TAU( I ) . . . TAUI24 ) ) 

Si) A - -<>99* 

VF- 1 

L WK.tDEPl. ) 

READ A.siO PRINT UATA 
FUNC TI ON RETURN 

** STEP 1 - INPUT COEFFICIENT FILE, IF NECESSARY ** 

** " INPUT PRCGRAM VARIABLES A \0 CALCULATE ** 

** NECESSARY ARRAYS 

S T £ P ( 2 ) 2 E R U * (CHI N»C FILE r XM I J , DE L f MPA I R$,PAlKS.. .PAIRS! 36),lCT«n 

PMOUE - t K 4 D J A N $ 

PROG = tDEPTHi 

2ERU. (CGfcFSw f OF£P| 1 ). . . Dh t P ( l 4 ) \ 

READ AND PRINT DATA 

WHENEVER COEFSW.fJF . D t TRANSFER TO bKO 
WHENEVER TFLAG.c'.O, TRANSPtK TO BR 1 
** READ IN A COEFFICIENT FILE *♦ 

WHENEVER CaiNl . t- -l - ANU.CGI N. E. 0 

PRINT FORMAT $H*0Ml COEFFICIENT !AP£ SPEC l F 1 EO« * % 

ERR ERROR. 

END OF CONDITIONAL 
** TAPE ACTION ** 

WHENCVER CHlN.E.CHlNl 
REWIND TAPE CUNIT 
OTHERWISE 

MOUNT* ICBlNtCUM T»$GUT*) 

CtWNl = C41N 
£ MO OF CONDITIONAL 
SKIP. (CFILE-l « 0 f CUN 1 T ) 

** READ FIRST RECORD OF COEFFICIENT FILE ** 

R BIN. (CUNI T , 1 A (2 ) t 3) 

KC = RCHECK • (0) 

IoHEnEVCR R C • is E • 0 

PRINT FORMAT KERR t RC f l 
TRANSFER TO ERR 
END UF CONDITIONAL 

** CHECK FOR ADDITIVE COEFFICIENTS ** 

N = QMS S * (JNCHAN 
WHENEVER I A * E . 0 

WHENEVER N.NE. (1A( 1)M A<2) I 

PRINT FORMAT *H • 0 1 NC UN S I S T LNCY BETWEEN COEFFICIENT FILE AND DATA F 

1 !IE*/H'0** COEF FILE ** NSS = * r 1 5 * S5 , H 1 NCHAN -SIS/H 1 ** DATA 

2 FILE ** NSS -• , I NCHAN = f , i b/ * $ * l A 11 ) , I A < 2 ) , QNS S , 

3 GNCHAN 
TRANSFER TO ERR 

END OF CONDITIONAL 
OTHER W l Sc 

PRINT COMMENT &ONON-AD01 T 1 VE COEFFICIENT FILE* 

TRANSFER TO ERR 
FND OF CONDI riONAL 
** ^EAO COEFFICIENT record ** 

K8IN. (CUN I T t ACOEF(N) ,NJ 
RC * RCHECK. 10) 

WHENEVER rc.ne.o 

PRINT FORMAT RERR,RC,2 
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flKO SETOM. < AC Of: Ft L. . .GNSSt 1 . . . UNCH AN ) 

W HtMVtK ClitFSrf.Nfi.O 

THROUGH C L OOP i , FUR I = 1 , 1 , I . 0 . GNS S 
( J = 1 , 1 , J.G.GNCHAN f /VCGfcFU.J) - -s)cEP(J) ) 
CL 00*1 C j-'v T I OOF 

L N 0 OP C J vOl TlC!W\l 


** CHECK IF iE* F-FUNCHUN MUST l3F; COM PUT FC ** 

WHENEVER TF L AG . L . 0 • A NO ♦ XM l U * L * 0 * TRANSFER TC bR l 
** COMPUTE F-FUNCHON ** 

WHCNE VER S U A • L • - $ 6 0 . 

PRINT COMMENT 60 SUN DEPRESSION ANCLE NCJ SPECIFIED FOR NEK 
1 DATA FILE* 

TRANSFER TO ERR 
F N U OF CrHOIUONAL 
wHE.NL VCR DtL . t ♦ 0 

WHENEVER OUANG.E.G* 

PRINT FORMAT .5H 1 DURABLE TO COMPUTE F-FUNCT I GNS « ** 

TRANSFFR TO ERR 
f. M 0 OF CONDITIONAL 
(ITHFPwl SL 

UOANG - DEL * PiSC 

OH A NO = -HXMC - GNAJ/CKP1 * ODANG 
END OF CONDITIONAL 
0=1./ COS. ( SUA*P1B0) 

X = OiVANG 

THROUGH LCOPltFGR I = 1 , 1 , I . G . C\SS 
F ( 2 ) = 1 » / i I . /COS. IX) + 0 ) 

LOOP l X = X + CJDANG 

PRINT RESULTS F 1 1 ) . . . F J 4N55> 

** CHECK PAIRS FOR CONSISTENCY ** 

HR I THROUGH L0QP2 r F0R I -1 * 1 * I ~G« NPA IRS 

WHENEVER PAIRS! I f i) *G . QNCHAN. OR. PA IRS! It/) .G. WNCHAN 

PRINT COMMENT tOCHANNE L-PA J R COMPONENT GREATER THAN NUMBER 
1 OF CHANNELS ON INPUT T APE t 
TRANSFER TO FRR 
END OF CONDITIONAL 

** COMPUTE G-HJNC T I ON FOR EACH PAIR ** 

L00P2 GU) * 1. / t A LPH A (PAIRS! I * 1 ) ) - ALPHA < PA 1*$ U * 2 I)) 

PRINT RESULTS Gil) 

** RADIANCE POUF. CALCULATIONS ** 

WHENEVER PMi)G£*E « *RAD! Aft* 

THROUGH LOCP 3 f F UR K= 1 , 1 , K. G . NPA I RS 
I = PAIRSCKtl) 

J = PAiRSIKtZ) 

o = iannm*Hin*TAiM i M / <rhuij>*hi j)*T4u<jn*EK[ iweki j) 

LOOP 3 A(K> = F.LQC. (C) 

PRINT RESULTS Q t A ( 1 ) 

BRMOUE = ID 

** reflectancc mode calculations ** 

OTHER W I SE 

** CHECK VALIDITY OF SUN SENSOR REGION ** 

WHENEVER NA 5»G* NI3 S 

PRINT FORMAT 1>H • UNAS GREATER THAN NBS ** NAS =*,[4tH» NRS =',!***♦ 


®MAL fk®tM 
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STEP13) 


L00P4 


LOOPS 
HR 2 


NAS,NBS 


F R A :'4 S K F ,< TO £ R R 

OR L- N E V C R NA S * G . ( UNA + I UNS S-l ) * CK M ) 

pRrvr format «h'o$um sensor region not 
transfer to err 


A V A ILAI3LE* *4 


EMU UF CONDITIONAL 

NASL = [U\ 3 - CNA } / yKi J + 

l = l ONA + l S j- i ) P ) 
WHENEVER NBS.G.I t MBS = 
Wf3Sl = INt>S - LNM/QKP + 
6 RHODE - Otl 
END UF CONDITIONAL 

BP1 = ( -\A - UNA) / uKP + 1 

URP - KP / gxP 
( J s l , l , I .G.NC t UFACTRIH = 
OF LAG = 2 


1 

l 

1 


i/51.0 f C=FACr«U+25 >*-i ./51.0) 


FUNCTION RETURN 


** STCP 3 - LINE OPERA 1 10^ S ** 


k = pstart + e 

LCTR a LCTR ♦ 1 
BP = BP1 

** CALCULATE A— VECTOR FRU*i SUN SENSOR REG l ON. ** 
** IF REFLECTANCE MU DC ** 

WHENEVER BRNODEt TRANSFER TO BR2 
SPREAD. (-1 t M A X < 1) .. .MAX (UNCHAN) ) 

THROUGH LflOPA 9 FOR 1 = N A S 1 t 1 , I . G. NbS 1 
VAL . ( DA TA ( K ) , ! ) 

THROUGH LOOP A v F Oft J = l vl tJ.G.LNCHAN 
WHENEVER Rh r ( I ) .G.MAX { I ) , MAX(I) = RbT(I) 
CONTINUE 

( I 5 i » l f l . G . 0 NC HAN# F w A X ( I ) = M A X ( I ) J 

EXP SM T . { F K A X # V S S # wNC H A N ) 

THROUGH LOOMS, FOR N - 1 t 1 f N.G.!\j PAIRS 
I * PAIRS! N# 1 ) 

J = PA IRS <N, ? I 

0 = ( RHOU )*RHGSS( l ) *VS$l I ) *TAU ( U *T AU< l I l / 
IftHOI J)*RHOSS( J1*VSS< J)*TAU( J)#rAU( J) ) 

A(N) = FLOG. (g> 

FUNCTION RETURN 


INTERNAL FUNCTION DEPl. 

FLOATING POINT ZU8),MEAN,$DFV 
I = UATUMiED! TCH ) 

WHENEVER I .o.MI .OR. I.L.LO 

ZERO. (DA TUM< l ) . . .LA HJM(NC) ) 

TRANSFER TO DBK I 
END UF CONDITIONAL 
SUEV = 0. 

KEAN = 0. 

THROUGH 1 LOOP, FOR M= i , 1 # M „ G. NPA I H S 
I = PAIRS! K# l) 

J = PAIRS (M»2) 

V2 ^ OATUMIJ) + ACOEF ( BP t J ) 

VI - DATUM! n + ACULF<0P,I) 

**VF 15 LEVEL UF FLUCTUATIONS IN DEEP WATER VOLTAGE** 
WHENEVER V2.L. VF.OR. Vl.L.VF 
Z(M) = 10.1 


DEBHNAS PAGE IB 
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1 LOOP 


I LOOP 1 


ILQ092 
DbR 1 


OTHrRwi Zk 
X = \ f? 

X - X / VI 

w - t LOG « ( A I 

l i M ) * F < H P > *• G ( N> * < a ( v ) + Q } 

END OF CUND 1 1 ICNAL 

WHENr VER / { 0 ) . L. 0- , / 4 >i > = -J 

WfAs = K.V'N + Z(N') 

•v £ A N = ^lAN / M'-VAtRS 

** CALCuL Arc STANDARD DEVI AT 1 ON mnH | Lr SCALING DATA 
THROUGH ILCOPl ,FiJ« M - 1 , i .F.G.NPAlftS, 

0 = li * I - MOAN 

snev = snt v * q * o 

L» A T U ^ < N ] = Z(M) * CM + 1 « 

whenever NPMRS.F.I 
CATUMI 3 ) - 0 
OTHER w [ SF 

SuE V - SORT- ( SUtV/INPAMS-l )) 

OATUM1NPAIRS + 2 ) = SDfcV * CM 
END OF CONDITIONAL 

DA F UM I. \PA I* 5+ i ) = MEAN * CM 
THROUGH I L GIJP2 ♦ F OR IMtUr.G.NC 
WHENEVER DATUM ( i) .G.H i f UATUM| I) = Ml 
CON I I NUT 
HP = *P f DBP 
FUNCTION RETURN 
END OF FUNCTION 


INTERNAL FUNCTION VAL- (ARRAY, XX) 

THROUGH LOUP. FOR i l =0 » L , I I . E . ONCHAN 

YY = ( XX* I ) + CNCH AN + I I 

WORD = VY / 4 

DISP = ( 5- i YY-WQRD**) ) *9 

LOOP RCTUI + U - (ARRAY! WORD ) ♦ R S . D l $ P ) . A . M 1 

FUNCTION RETURN 
END OF FUNCTION 


INTERNAL FUNCTION EXPShT. UM«GUT'NN> 

FLOATING POINT I W , OUT f l l f l \ 

WHENEVER LCTw.G.lO 
TRANSFER TO ELCGP 
OR WHENEVER LCTR.F-10 
Ll - Oml 
LI = 0-9 

OR WHENEVER LCTR.G-l 
LL = U/LCTR 
1 i = i- - ll 
OTHERWISE 
Ll a 0- 
IV = l. 

END UF CONDITIONAL 

FLOOR (I = UUt-G-NNt OUT (I) * UU Ul}*ZZ ♦ INU)*ZL ) 

FUNCTION RETURN 
FNO OF FUNCTION 


END OF FUNCTION 


OMGWABPf®, 


Of 
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Appendix B 

COMPUTER PROGRAM FOR OPTIMUM DECISION- BOUNDARY METHOD 

A computer program also was written to implement the method of depth calculation dis- 
cussed in Section 3.4. This program can use any number of input channels and writes out the 
depth values at each point on a single-channel output tape. The input data include the signal at 
the known depth, which may be zero (using the highest signal observed at the water's edge). 

The input variables are summarized, in Table B-l, and the computer coding is listed in 
Table B-2. 
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TABLE B-l. VARIABLES FOR DEPTH2 PROGRAM 


Card 1: 



Input Variables 

Mode 

; - Description 

NV . 

Integer , , , ; 

No. of channels on input tape 

NCHAN 

Integer 

No. of channels used in the depth 
calculation (NCHAN < NV) 

CHAN (I) 

Integer 

Channel number used (1=1. . . N CHAN) 

Card 2 

(process variables): 



same as for DEPT HI 

. ■ . 


Card 3: 



ALPHA(I) 

Floating Point 

Water -attenuation coefficient in CHAN(I) 

DEEP(I) 

Integer 

Deep-water signal in CHAN(I) 

VREF(I) 

Integer. 

Reference signal in CHAN(I) 

ZREF 

Floating Point 

Water depth where reference signal is 
measured 

SDA 

Floating Point 

Solar -z enith angle 

EDITCH 

Integer 

Same as in DEPTH1 

LO 

Integer 

Same as in DEPTH1 

HI 

Integer 

Same as in DEPTH1 

• ' • ' * 


Notice that this program doe s : not require a knowledge of the bottom reflectance, sensitivity 
constants, solar irradiance, or atmospheric transmissivity. All the necessary information is 
contained in the reference signal VREF(I). ; 

Scan -angle variations have been neglected, eliminating the input variables DEL and XMID. 

If there are variations in bottom reflectance in the scene, the depth values will be in error 
by an additive factor for the areas where the bottom reflectance differs from that where 
VREF(I) was measured. In some cases, however, it might be possible to choose a set of 
channels to minimize or eliminate this error. 
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TABLE B-2. LISTING OF DEPTH2 PROGRAM 


EXTERNAL FUNCTION DLPTH2 « 

NORMAL NODE IS INTEGER 

references on 

ERASABLE SKP2(226), OUNIT, OREEL, CFILE, OLINE, THRITE, OCHAN, 

1 Oft Ml) S', PROG, KODE1, MCJUE2 , . U J I T , CAL I ME, COLINE, NS A, NSB, 

2 KS. NA, NCi , KP, lUlll, TFLAG, IPACK, RESERVIA), OFACTR(A9), 

3 CT1TL2II9) , CTITLU19), ELI SI ( 19 ) , CS PARE (AG ) , OFLAG, 

A OOANG, ObANG, CKP, UNA, QNSS, ONCHAN, QMOUc, uRECA , QRECD, 

5 LFILF, BREEL , CNWROS 

ERASABLE OATAIA23), ITEST , OSTART, OSTAP.T, NV , NX, NC, L, IP, 

1 IGP, NOP, NEXT, STAR1, REACT, REaUL, REaQP, LABEL, 

2 EXTRA (20) , DATUM ( 2A ) , IC0UEI2A), ICHANI2A), COEF(l) 


DIMENSION CHAM2A) ,VREF(2A) ,0EEP(2A) 

FLOATING POINT ALPH A ! 2 A ) , QFACTR , SDA, X, P 180 , C5 l 
FLOATING POINT l , ZR EF , ZM AX ,C , V 
FLOATING POI NT SORT. , COS. , ELOG. 

VECTOR VALUES Cbl=50.0 
VECTOR VALUES PI R0=0.017C5329 
WHENEVER NEXT. G. 2, FUNCTION RETURN 
TRANSFER TO STEP ( NEXT ) 

STEP ( 1 ) NC-1 

LINK. (0F.P2.) 

READ AND PRINT DATA N V , NCHAN.CHAN U ) . . . CHA.NI NCHAN I 
FUNCTION RETURN 

STEPC2) ZERO. (DEEP Cl).,. DEEP (2A) > 

ZERO. (EDI TCH, LO.HI ) 

ZKEF=0. 

READ AND PRINT DATA ALPHA ( I) ,D£EPl I) ,V REF ( 1 ), ZREF, SDA, EO ITCH, LO.HI 
C = 0. 

THROUGH DRL1, FOR I = l , 1, I . G.NCH AN 
0RL1 C=C+ALPHA( I) * ALPHA) I) 

C=C* ( l. + l./COS. ISDA*P180) ) 

ZKAX=ZREF 

THROUGH 0RL2 , FOR 1 = 1, 1, 1. G. NCHAN 
V= VREF ( I ) -DEEP ( I ) 

DRL2 ZMAX=ZMAX+ALPHA( I )*ELOG. IV)/C 

PRINT RESULTS ZMAX 
■iFACTSI L I =1. /C51 
OFACTR (261=0. 

0FLAG=2 

FUNCTION RETURN 

INTERNAL FUNCTION 0EP2. 

I = DATUM (ED ITCH) 

WHENEVER I.G.HI.OR. I.L.LO 
DATUM (11=511 
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DRL 3 
DRL4 


TRANSFER TU URL4 
END OF CONDI T I GNAL 
L - ZHAX 

THROUGH DHL 3 » FOR l =1 » 1 , I . G. NCH AN 
V=DATUMICHA‘HI ))-OEEP< 1) 

WHENEVER V.LE.O. , TRANSFER TU UKL3 
X = ELOr,. (V) 

Z ~ Z - ALPHA < l I*X/C 
CONUNUt 

WHENEVER Z.L.O., Z=G. 
l)A TUMI 1 I = Z*C5 1+0 . 5 

WHENEVER DATUM (1 1 .G.500 ,DATUrtU. ) = 500 

CONTINUE 

FUNCTION RETURN 

ENO OF FUNCTION 

END OF FUNCTION 


PAGE W 
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